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ABSTRACT 

Co-addition of blank-field FIRST data at the location of over 8000 QSOs in the 2QZ 
survey has yielded statistical detections of radio quiet QSOs with median flux levels of 
20 — 40/iJy. We show that the total radio flux of radio-quiet QSOs in the 2QZ is consis- 
tent with a smooth extrapolation of the 2QZ radio-loud QSO number-flux distribution 
based on the slope of the relation flattening near the FIRST flux limit. However, we 
are unable to distinguish a smooth extrapolation of the luminosity function to faint 
levels from a bimodal luminosity function with a break close to or below the FIRST 
radio detection limit. We also demonstrate that the redshift dependence of the median 
radio-to-optical spectral index chro for these radio quiet QSOs is consistent with that 
obtained for individual radio-loud 2QZ QSOs detected by FIRST. 
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1 BACKGROUND 

The debate over whether the observed properties of radio- 
loud and radio-quiet QSOs are consistent with a single pop- 
ulation of objects, albeit with a broad range in radio proper- 
ties, or two distinct populations is more than thirty years old 
(see e.g. Schmidt 1970). Evidence for bi-modality in the ra- 
dio properties of QSOs claimed by a number of authors (see 
e.g. Kellerman et al. 1989, Miller, Peacock & Mead 1990) 
has been challenged more recently using data obtained from 
the deeper, more extensive radio and optically-selected cat- 
alogues (White et al. 2000, Hewett, Foltz & Chaffee 2001) 
that are increasingly becoming available. However, recent 
claims for bi-modality have also been made based on QSO 
samples drawn from the extensive Sloan Digital Sky Survey 
(Ivcsic et al. 2002, although see Ivesic et al. 2004 for a re- 
assessment of this result) In contrast, a statistical analysis of 
the 2dF QSO redshift survey (2QZ, see Croom et al. 2004) 
and LBQS optically-selected survey and FIRST radio sur- 
veys has yielded no evidence for any gap in the distribution 
between radio-loud and radio-quiet populations (Cirasuolo 
et al. 2003). 

Due to the instrincally low radio flux of radio-quiet QSOs 
<C 100/iJy, most previous analyses have been based on 
statistical arguments applied to upper flux limits inferred 
from non-dectections, as opposed to detection, of radio quiet 
QSOs. In this paper, we adopt a different approach by di- 
rectly trying to obtain a detection for an 'ensemble' aver- 
age of radio-quiet QSOs in the 2dF QSO survey by stack- 



ing radio data from the FIRST survey. The combination of 
the 2QZ sample size (~ 10000 QSOs in the area of overlap 
with FIRST) and the moderate depth of the FIRSF sur- 
vey (0.15mJy rms) suggest it should be possible to achieve 
reliable detections (> 5cr) for a number of 1000-2000 QSO 
ensembles at a level of a few tens of fijy. Indeed, this method 
has already been successfully pioneered on the Sloan QSO 
sample by Glickman et al. (2004). 

We describe the data and method used in Section 2 and 
present our analysis of the radio-quiet QSO flux distribu- 
tion derived directly from the coadded images in Section 3. 
'Standard' cosmological parameters: f2 m = 0.3, £Ia = 0.7 
and Ho — 70kms~ 1 Mpc~ 1 are used throughout this paper. 



2 METHOD 

We used the 2QZ to obtain accurate (0.2-0.3 arcsec) posi- 
tions for 8741 optically-selected QSOs within the 75-degree 
x 5-degree wide 2QZ equatorial strip (9 h 40 m < RA < 
14 h 40) that overlaps with the FIRST survey (White et al. 
2001). To ensure maximum data integrity, we used only 
QSOs with ID and redshift quality 1 in the 2QZ, and re- 
stricted our analysis to QSOs in the redshift range z < 2.3 
where the survey completeness is over 50 per cent (Croom 
et al. 2004). We did not exclude the lowest redshift QSOs 
(z < 0.25) from the analysis where the completeness is also 
less than 50 per cent. However these objects only comprise 3 
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Table 1. Radio properties of composite QSOs 
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Figure 1. Images from co-added redshift bins. Increasing redshift 
bins run from left to right and top to bottom. The image in the 
lower right part of the figure is the sum of all redshift bins. 




per cent of the QSOs in the lowest redshift bin under study 
here, and we have verified that their inclusion/exclusion 
makes no difference to the results obtained below. 
Out of this sample of QSOs, 228 are listed in the 2QZ pub- 
lic release catalogue as having NVSS counterparts. A further 
134 were also found to have FIRST radio counterparts; de- 
fined as a source listed in the FIRST catalogue (i.e. > 5a 
detection) within 5 arcsec of the optical position. The FIRST 
radio fluxes for all these 362 QSOs, hereinafter referred 
to as radio-loud, lay in the range 0.65mJy< /(1.4GHz) < 
3092 mjy. All 2QZ objects identified as a radio-loud quasar 
were consistent with a point source radio detection. 
For the remaining 8379 QSOs we extracted 30 arcsec x 
30 arcsec images centred on the optical position from the 
FIRST catalgoue. Each 'cut out' was visually inspected, and 
70 were rejected from further analysis based on either poor 
data (sidelobes etc.) or the presence of a bright source in 
the edge of the field. The remaining 8309 QSOs, hereinafter 
radio quiet QSOs, were then grouped into five redshift bins. 
Each bin was designed to have an equal number of QSOs, 
so that the rms noise in each stacked image would be simi- 
lar. The redshift bins and the numbers of QSOs are listed in 
Table In practice, there is a small variation in the num- 
bers of QSOs between each bin, due to an iteration in the 
removal of maps with poor data. 

The radio images of the QSOs in each redshift bin were 
stacked using the IRAF imcombine routine using the me- 
dian estimator (see Fig. 0. As a control sample, we also 
stacked the similar FIRST 'cut outs' around the positions 
of 4830 galactic subdwarfs and 885 white dwarfs in the 2QZ. 
Galactic subdwarfs are not expected to exhibit radio emis- 
sion, even at the few /iJy rms noise levels attained with a 
stacks of over 4000 0.15mJy-rms images (see Fig|2J. 
The imfit routine in MIRIAD was then used to determine 
the integrated flux of each source detected in the stacked 



Figure 2. Stacked images from both white dwarfs (left image) 
and galactic subdwarfs (right), from FIRST cutouts at positions 
from 2QZ. 




images. Images were detected at greater than the 5<r level 
in each co-added QSO image, but no source detection was 
made in either the galactic subdwarf or white dwarf control 
stacks. Images were modelled as point sources, and a gaus- 
sian fit was used with a 4-pixel (7.2-arcsec) diameter aper- 
ture. This aperture is consistent with the 5-arcsec resolution 
of the FIRST survey and with the Gaussian 5-arcsec FWHM 
aperture used to derive the FIRST point-source fluxes used 
for the radio-loud objects in this analysis. The flux estimate 
was robust against the choice of the precise centring (±0.5 
pix) of the aperture. Furthermore, increasing the radius of 
the aperture produced no significant change in the total de- 
rived flux. 

The median QSOs in the redshift bins were found to range 
in flux from /(1.4GHz) mo( j = 23/iJy to 42/xJy. This point- 
source flux corresponds to angular scales of between 20 kpc 
(z = 0.65) and 45 kpc (z = 2.05) for the median redshift 
bins used in this analysis. Given that, the vast majority of 
quasars is this survey will be compact flat-spectrum sources 
(see below) this 'point source' flux should represent an unbi- 
ased estimate of the total radio flux, independent of redshift. 
We note that derived fluxes for these stacked sources are 
similar to those obtained in the similar analysis made by 
Glickman et al. (2004) on the stacked radio-quiet QSOs in 
the Sloan QSO survey. 

3 ANALYSIS 

3.1 Number counts 

We first calculated the observed radio number-flux rela- 
tion, n(s), for both the radio- loud and radio-quiet QSOs 
in the 2QZ to determine whether there is evidence for any 
bi-modality or discontinuity in the flux distribution which 
might be indicative of a two-population model. We used the 
catalogued FIRST 20cm fluxes for the individual radio-loud 
2QZ QSOs. The resultant n(s) is shown in Fig [3] 
We wish to test whether the observed total flux in the 
stacked images is inconsistent with a smooth extrapolation 
of the observed radio-loud n(s). We know that the total 
flux of the radio-quiet population can not exceed the to- 
tal flux, /(1.4GHz) t ot obtained by summing the median 
fluxes detected in each of stacked images over all redshift 
bins, scaled by the number of QSOs in each stack i.e. 
/(1.4GHz) tot = Eti ^QSoW x /(1.4GHz) mcd (i). 
On the basis that the radio n(s) may be extended below the 
FIRST completeness limit, we can obtain a prediction for 
the total QSO radio flux based on a smooth extrapolation 
of a fit to the radio n(s) above the FIRST completeness 
limit. We adopt a completeness limit of 5<r (7.5 mjy) for 
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the FIRST survey, and fit the the n(s) between the 7.5 mjy 
and 20mJy. Compared to higher fluxes, the slope of n(s) 
relation for radio-loud QSOs 'flattens-off ' at these flux levels. 
Fitting of the n(s) relation suggest the radio counts may 
be extrapolated below 7.5 mjy by a fit ranging from a flat 
slope, n(s) — constant, to shallow slope n(s) oc s -0 ' 30 , with 
a best-fit n(s) oc s" 015 (see Fig. EJ. 

Table [5] presents the comparison of the integrated flux as a 
function of depth for these fits between these models and 
stacked observations. The errors on the integrated flux in 
both radio-loud domain and radio-quiet domain are approx- 
imately 10 per cent. Both model extrapolations converge 
rapidly at fluxes less than O.fmJy, with integrated flux pre- 
dictions for slopes of 0.15 and 0.30 bracketing the observed 
integrated flux from the 2QZ radio-quiet QSOs. We thus 
conclude that the flux distribution of radio-quiet QSOs in 
the 2QZ is not inconsistent with a smooth extrapolation of 
the 2QZ radio-loud QSO number-flux distribution, on the 
basis that the slope of the relation flattens at faint flux lev- 
els. 

This is a necessary but not sufficient condition for a sin- 
gle population model; more detailed modelling of the source 
counts based on fitted radio-loud and radio-quiet LFs will be 
presented elsewhere (Miller et al. in preparation). We note 
that the cosmological luminosity evolution of QSOs in both 
optical (Croom et al. 2004) and radio (Dunlop & Peacock 
2004) tends to cancel out the effects of cosmological dimming 
over a significant range in redshift (0.5 < z < 2). This results 
in luminosity function features such as the break, L*(z), ex- 
hibiting the same observed flux over these redshifts. Thus a 
simple bimodal luminosity function model with a redshift- 
dependent break between radio-loud and radio-quiet QSOs 
that maintains an observed flux close to or below the FIRST 
radio detection limit over the range in redshifts sampled here 
could reproduce the current observations. 
We also note that the two-population model favoured by 
Miller et al. (1990) was based on a carefully-selected sample 
of QSOs at 1.8 < 2 < 2.5. We therefore looked to see if there 
was any evidence in any of our redshift ranges for the inte- 
grated radio-quiet flux below the FIRST completeness limit 
departing from a smooth extrapolation of the radio-loud 
n(s). Due to reduced numbers of objects, the faint slope was 
much more difficult to constrain in these narrower redshift 
bins and we retained the values of 0.0 and 0.3 derived for the 
full n(s) above. The results are reported in Table based 
on an integration from the FIRST 5a limit down to 10/iJy 
. As an illustration the n(s) relation for the 1.85 < z < 2.30 
range is plotted in Fig. Although hampered by larger 
measurement uncertainties, (~ 20 per cent), we found the 
total integrated observed flux below the FIRST complete- 
ness limits (0.35 — 0.48 mjy deg -2 ) at these redshifts were 
consistent with our extrapolation of the radio n(s). In the 
lowest redshift bin, the total observed flux marginally ex- 
ceeds the predicted flux with the fitted 0.3 slope, but only 
at the lcr level. 



3.2 Spectral Index 

We also investigated the dependency on redshift of the 
radio-to-optical luminosity ratio for the radio-loud and 
radio-quiet QSO populations. We characterised the radio- 
to-optical luminosity ratio by a IO , the notional spec- 
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Figure 3. Number-flux relation, n(s) for the 2QZ radio-loud 
(filled squares). Fits to the n(s) used in the text are indicated by 
the dashed lines (slopes of 0.0, 0.15 and 0.30 respectively). The 
5(T completeness limit and 20 mjy fitting limit for the FIRST 
radio survey are shown by the dotted lines. Also shown in the 
corresponding n[s) relation for the range 1.89 < z < 2.30 (open 
squares). 
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Figure 4. Radio to optical spectral index as a function of redshift 
for both radio- loud (crosses) and radio-quiet QSOs (filled circles). 



tral index between the radio and optical bands, a ro = 
— (logZ/440nm/ii.4GHz)/(log^440nm/^i.4GHz) where L refers 
to the luminosities in both optical (440nm) and radio 
(1.4GHz) bands. Luminosities were calculated using a spec- 
tral index of —0.5 in both the optical and radio bands. The 
assumption that the mean radio and optical spectral indices 
are independent of the radio 'loudness' is implicit in our cal- 
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Table 2. Integrated Flux density for observed QSOs and extrapolated n(s) 



Integrated Flux Density 
Flux Limit n(s) = const n(s) oc s -013 °c s~ 0,2 Observed 

(mjy) (mjy deg -2 ) (mjy deg -2 ) (mjy deg -2 ) (mjy deg -2 ) 



5.0 


0.50 
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0.76 


0.64 


3.2 


0.82 


1.03 


1.31 
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2.0 


1.02 


1.32 


1.71 


1.13 


1.2 


1.14 


1.50 
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0.8 


1.22 


1.63 


2.20 


1.27 


0.08 


1.35 


1.87 


2.64 


1.27 


0.05 


1.35 


1.88 


2.67 


1.52 


0.03 


1.36 


1.89 


2.70 


2.00 


0.02 


1.36 


1.90 


2.71 


2.10 


0.01 


1.36 


1.90 
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2.10 



Table 3. Observed and predicted integrated flux densities at f(1.4GHz)> 0.01 mjy for different redshift bins 



Integrated Flux Density 
z range n(s) = const n(s) oc s — 3 Observed 
(mjy deg -2 ) (mjy deg -2 ) (mjy deg -2 ) 



0.12 < 2 < 0.85 


0.24 


0.40 


0.48 


0.85 < 2 < 1.23 


0.30 


0.50 


0.40 


1.23 < 2 < 1.55 


0.32 


0.53 


0.38 


1.55 < 2 < 1.89 


0.32 


0.53 


0.34 


1.89 < 2 < 2.30 


0.26 


0.43 


0.38 



culations. Studies of composite optical spectra for radio-loud 
and radio-quiet quasars confirm that the optical spectral in- 
dex in the range 1200 A < A4400A is very similar for both 
populations (Cristiani & Vio 1990, Zheng et al. 1997). Al- 
though there is no information of the radio spectral index of 
radio-quiet QSOs, at the redshifts under consideration here, 
evolutionary models of the radio source population suggest 
that the vast majority of the objects will be flat radio spec- 
trum QSOs largely independent of redshift (see e.g. Jackson 
& Wall 1999); steep spectrum sources only dominating at 
low redshifts z < 0.5. 

In Fig 2] the derived a IO for radio- loud and the composite 
radio-quiet QSOs are shown. Although the absolute redshift 
trend in a IO is dominated by the choice of spectral index, the 
relative trend between the radio-loud and radio-quiet QSOs 
(which is independent of radio or optical spectral index) is 
similar. This is confirmed in Figure |S] where we plot the dif- 
ference in Q ro between the radio-loud and radio-quiet QSOs 
as a function of redshift. Error bars are determined from the 
errors on the flux determinations in both optical and radio. 
We can express the relation thus: 

Lrq(z) _ Lrl(z) 
Prq{z) ~ 'Prl(z) 

where the subscripts refer to the radio-loud (RL) and radio- 
quiet (RQ) populations respectively and, for clarity, we de- 
note the optical luminosity (Z/440nm) and radio luminosity 
(ii.4GHz) as L and P respectively. 

Over the redshift range of interest (z < 2.3) the observed 
radio flux evolution of the radio-loud population and the the 



optical flux evolution of the radio-quiet populations can be 
parameterised by luminosity evolution model as follows: 

Lrq(z) = Z / *i0 felz+fc2z2 

and 

Prl{z) = p*i0 ft 'i z+fc ^ 2 

where the optical evolution fit to the radio-quiet population 
yields ki = 1.39 and k<z = —0.29 (Croom et al. 2004), similar 
to the radio evolution fit to radio-loud QSOs fei = 1.18 and 
ki — —0.28 (Dunlop & Peacock 1990). By combining these 
relations and expressing the functional forms for the optical 
and radio evolution as f op t(z) and / ra d(z) respectively, we 
obtain: 

Prq(z).L rl (z) = const./ ra d(«)./opt(«) 

We can use the observation by Padovani (1993) that the op- 
tical evolution of radio-loud and radio-quiet QSOs are iden- 
tical (Lrl(z) = Lrq(z)) to yield: 

Prq(z) = const. / tw i(z) 

i.e. that the radio evolution of radio quiet QSOs is identi- 
cal to the radio evolution of radio loud QSOs. As indicated 
above, this is predicated on the assumption that the mean 
optical and radio spectral indices of QSOs are independent 
of their radio-loudness. However, if this were not the case, 
the observation result would itself rely on some fine tuning 
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Figure 5. The difference between [a ro ]flL and [ct ro ]RQ, as a 
function of redshift. 

i.e. the relative difference in spectral indices between radio- 
loud and radio-quiet QSOs would have to be equal and in 
the opposite sense to the difference in evolution between the 
radio-quiet and radio-loud populations. 
This is consistent with the results of Glickman et al. (2004) 
who found that the radio-to-optical luminosity ratio of the 
stacked radio-quiet QSOs in the SDSS sample exhibited the 
same dependence on redshift as the radio-loud SDSS QSOs. 



4 CONCLUSIONS 

We have demonstrated that the co-addition of blank fields 
in the FIRST survey provides a robust detection of an en- 
semble average of radio-quiet QSO flux distribution below 
the FIRST threshold, and that the flux properties are con- 
sistent with a smooth extrapolation of the radio-loud QSO 
flux distribution detected above this threshold. 
It would be straightforward to apply this technique to other 
extensive deep optical surveys, (2dFGRS, SDSS galaxies) 
comprising tens to hundreds of thousands of objects, thus 
potentially moving into the sub-/iJy detection limit. With 
the advent of next generation wide-field radio facilities such 
as the Square Kilometre Array, the nanojansky detection 
limits proposed for surveys could be similarly transformed 
into sub-nanoJansky limits using these stacking techniques, 
providing the confusion limit is not reached and suitable 
object catalogues exist at deep levels in other passbands. 
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